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Abstract 

Experimental data have provided stringent constraints on neutrino mixing parameters. In the standard 
parameterization the mixing angle #23 is close to 7r/4. There are also evidences show that the CP violating 
phase is close to — 7t/2. We study neutrino mass matrix reconstructed using this information and find 
several interesting properties. We show that a theoretical model based on the A 4 symmetry naturally 
predicts <5 = —7r/2 and $23 = vr/4 when the Yukawa couplings and scalar vacuum expectation values are 
real reaching a y — r exchange and CP conjugate symmetry limit. In this case CP violation solely comes 
from the complex group theoretical Clebsh-Gordan coefficients. The model also predicts | V^ 2 1 = l/\/3 
consistent with data. With complex Yukawa couplings the values for 5 and $23 can be significantly deviate 
away from the symmetry values —7r/2 and vr/4, respectively. But \ V e 2 \ = l/\/3 is not altered. This matrix 
is an excellent lowest order approximation for theoretical model buildings of neutrino mass matrix. 
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Tremendous experimental progresses have been made in obtaining information about the neu¬ 
trino mixing parameters. The mixing angles in the Pontecorvo-Maki-Nagakawa-Sakata Vpmns 
matrix[l] are not always small [20] • In the standard parameterization[2, 5| for three neutrino mix¬ 
ing commonly used|$, 0], the mixing angle 0 23 is close to 7 t/ 4 , d 12 is large, d 13 is relatively small but 
away from zero, and also S 12 C 13 is close to l/\/3. Since the mixing angle d 13 is non-zero, the famous 
tri-bimaximal mixingis ruled out. There are now evidences show that the CP violating phase 
6 is close to — tt/2. This also implies that the tri-bimaximal mixing is in trouble since it predicts 
5 = 0. The phase 5 is sometimes referred as Dirac phase which shows up in neutrino oscillations. 
If neutrinos are Majorana particles, there are also new CP violating Majorana phases cq. There 
are many discussions about implications for data available emphasizing the particular values for 
|h| = 7 t /2 and 6*23 = vr/4[8-10]. One of the commonly mentioned property for this type of mixing 
is the so called maximal CP violation because |5| is 7 t/ 2 . This is, strictly speaking, an incorrect 
statement because that the value of the Dirac phase is parametrization dependent. For example, 
even the absolute value of the Dirac phase is 7 t /2 in the standard parametrization, in the original 
Kobayashi-Maskawa parametrization for quarks [ll| it is not 7 t /2 anymore. However, the special 
values for some of the mixing angles and the Dirac phase can still provide important information 
about neutrino mass matrix and can guide theoretical model buildings to search for the underlying 
theory. 

To this end, let us reconstruct the neutrino mass matrix assuming that neutrinos are Majorana 
particles with 8 = —tt/2 and d 2 3 = vr/4. In the basis where the charged lepton mass matrix is 
already diaganolized, the neutrino mass matrix defined by the term giving neutrino masses in the 
Lagrangian (1/2has the following form 


m u — VpMNsVivVpMNS , 


( 1 ) 


where m u = diag(mi,m 2 ,m 3 ) with m, = \rrii\exp(iai). Here we have put Majorana phase informa¬ 
tion in the neutrino masses. The standard form for Vpmns is given by 


C 12 C 13 

VpMNS = I — s 12 c 23 — c 12 s 23 s 13 e 


iS 

iS 


S 12 C 13 
C 12 C 23 — S 12 S 23 S 13 e 


i5 

.A 6 


’ 12 0 23 


'12 u 23 ms 1 


'12 °23 


’12 °23°13 1 


S 13 e " 
s 23 C 13 


'23 m3 


( 2 ) 


where Cij and s 


following form [ 8 , 10 


are cos dij and sin 0 l3 , respectively. With 8 = —7 t/2 and 6*23 = 7 t/4 , m v has the 


a c + i/3 —(c — i/3 ) 

m, y =| c + i/3 d + iy b 

— (c — i/3) b d — i 7 


(3) 


where 


a = rmc 2 12 c 2 13 + m 2 s 2 12 c 2 13 - m 3 s 2 13 , 


c =-+r( m l ~ "D)Si 2 Ci 2 Ci 3 , 

v 2 


b— ^ + c 12 s 13 ) + m 2 (c 12 + s 12 s 13 ) m 3 c 13 ) , 


d—^{mi(s 12 c 12 s 13 ) + m 2 (c 12 s 12 s 13 ) + m 3 c 13 ) , 


P = jyf Sl3Cl3 ( m l C 12 + m 2 S 12 + m 3) , 7 = -("H - "l 2 )Sl 2 Cl 2 Si 3 . 


(4) 


Note that in the most general case, because non-zero Majorana phases, the parameters a, b, c, <3, 
/3 and 7 are all complex. 
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The above matrix has a high level regularity pattern implying some underlying symmetry may be 
at work to produce it. Searching an underlying theory guided by symmetry principle may achieve 
this. Before doing this, however, it is worthwhile to understand more about the mass matrix in 
eq. (15|) . An immediate question one may ask is that if, in general, the neutrino mass matrix in 
eq. (13T) always predicts 8 = —7t/2 and d 2 3 = vr/4. The answer is negative. If 8 = 7t/2 and #23 = 7t/4, 
the neutrino mass matrix is given in a similar form as that in eq. (j5|) . but f3 and 7 need to be 
multiplied by a sign. Therefore without further information given, a general mass matrix in 
the form given by eq. ((3j) can give 5 = ±7t/2 and 6*23 = 7t/4. Whether they predict +7t/2 or —n/2, 
additional information need to be provided. Moreover, If neutrinos have Majorana phases, the 
general form does not imply that 8 and 623 must take ±7t/2 and 7t/4, respectively, neither. This 
can be understood by studying the following quantity 

m u ml = VpMNsrri v mlV ] PMNS . ( 5 ) 

The general form for neutrino mass in eq. (fill) will give the “12” and “13” entries A 12 j i 3 of m v m\, as 
Ai 2 + Ai 3 = — i2(a[3* + cq* — /3d* — /3b*) 

= -(|mi | 2 - |m 2 | 2 )si2Ci2Ci3(c23 - S 23 ) 


A12 — A 13 — 


13 


-(|mi| c 12 + \m 2 \ s 12 - \m 3 \ )si3Ci 3 (c 2 3 + s 2 3)e 
= 2 (ac* + cd* — cb* + /3q*) 

= -(|mi| 2 - |m 2 | 2 )si2Ci2Ci3(c23 + S 23 ) 

+ (|mi| 2 c 2 2 + |m 2 | 2 s 2 2 - |m 3 | 2 )si3Ci 3 (c23 - S23)e 


-iS 


-iS 


( 6 ) 


If the parameters in the set P : {a, b, c, /3, 7 }, are complex, the above equations can find solutions 
for other values of 623 and 8 . Therefore the general neutrino mass matrix form does not imply that 
8 and 0 2 3 must be ±7t/2 and 7t/4. If, however, the parameters in the set P are all real, as long as 
sin 8 ^ 0, one must have s 2 3 = c 2 3 and 8 = ±7t/2 as can be seen from the above two equations. 
From eq. TO) and eq. fRJl) . one also finds that all eigen-masses m; are real (the Majorana phases are 
zero or n). In this case the neutrino mass matrix can be rewritten as 


m v = 



(7) 


with A = a, B = b,C = c + i/3, and D = d — iq. The most general m u can be written as 


m,, = 


oipi 

0 

0 


0 

JP2 

0 



0 

JP 2 

0 


0 

0 

7p3 


( 8 ) 


where the phases pi are arbitrary. 

All neutrino mass matrices which can be written in the above form, will predict 8 = ± 7 t/ 2 , 
623 = vr/4 and all the eigen-masses are real. One can choose some particular values for p^ to obtain 
forms of ra„ for convenience of analysis. For example the “-” sign for the “13” and “31” entries can 
be removed by choosing p\ = p 2 = 0 and P 3 = n, the resultant matrix can be written in a more 
familiar forms 


ACC* 

m v = I C D* B 1 , 
C* B D 


(9) 
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where B = —B. 

The simplicity of the above mass matrix may serve as a good starting point to understand the 
possible underlying theory. If this has something to do with reality, one should not stay at the pure 
phenomenological level for analysis, but go further to study whether there are theoretical models 
which can obtained such a neutrino mass matrix in some consistently way. Several attempts for 
model buildings have been madejs], [9|. It has been shown in ref. [7f by Grimus and Lavoura that 
the above form of mass matrix is symmetric under a transformation of e —> e, p — r exchange with 
a CP conjugation. We will refer this as the Grimus-Lavoura symmetry (GLS). In this work we 
start with a simple model proposed earlier based on A4 symmetry fl3l] to realize the tri-bimaximal 
neutrino mixing, and then modify it to allow a non-zero 0 13 to find the conditions for having the 
GLS limit for neutrino mass matrix with <5 = —7 t/ 2 and #23 = vr/4 and how modifications may 
occur by explicit model studies. This model has an added bonus thatfl2l. 0 Sl2Cl3 — V e 2 — 1/\/3. 
There is also an interesting feature in this model that CP violation can be solely from complexity of 
relevant Clebsh-Gordan (C-G) coefficients in the GLS limit. We will refer this property as intrinsic 
CP violation. 

In this model A4 is serving as a family symmetry [l3| . The Higgs sector is enlarged to have three 
Higgs fields, $ = ($ 1 , $ 2 , $ 3 ) (SM doublet), 0 (SM doublet) and x — (xi, X 2 , Xs) (SM singlet). 
Under the A4, $ and x both transform as 3, and 0 as 1. Three right-handed SM singlet neutrinos 
vr = (isr, u%, iSft) are introduced allowing seesaw mechanism to be in effective. The standard 
left-handed leptons U = (l\, l\)i and standard right-handed charged leptons {Ir,Ir,Ir), and 
vr transform as a 3 , (1, V'. P) and 3, respectively. We refer the readers for more details on A 4 
group properties to Refs. (^. Il3l. lij . The Lagrangian responsible for the lepton mass matrix is 


L = A e (Z“ L $) 1 Z^ + A M (Z L $) r Z| + A T (Z L $) 1 «Z| + ^.C'. 

+ A„(/ l z/ r )i0 + m(v R v%) 1 + A x (P r z/^) 3 x, (10) 

where 

(Wik =(h*i + 7L*2 + il* 3 )i 1 R , 

(hfyl'l 1 * = + ^ 3 )ll , (11) 

(h$)l"lR = (1 1 l $i+u J 2 P l $2+cuP l $ 3 )1 3 r . 


Here cu = exp(i27i/3) and c o 2 = expiiAit/3) are the C-G coefficients of the A 4 group products. 

If the vev structure is of the form < $ 1 , 2,3 >= < Xi ,3 >= 0, < y 2 >= v x , and < 0 >= v^, 

one would obtain the charged lepton mass term as 

/ V3\ e v$ 0 0 \ / l l R \ /1 1 1 \ 

(ll ll A ) Ut 0 V3\»v* 0 h| , Ui = -J= 1 a; a; 2 . (12) 

\0 0 \/3A r n$ / \ l R J v d y 1 UJ 2 LO J 

From the above, we can identify the charged lepton mass to be m* = \/3AThe neutrino mass 
matrix has the seesaw form with 


M 


(0 m d \ 

V Ml Mr ) 


M R = 



0 rn x \ 

m 0 J , 

0 ml 


(13) 
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where Mr> = Diag(l, 1, l)X v v^, and m x = X x v x . From this one obtains the light neutrino mass 
matrix M u of the form given by 

( w 0 x \ 

0 y 0 , (14) 

x 0 z J 


where w = z = — (A u v^) 2 m/(m 2 — m ^), x = (X u v c f > ) 2 m x /(m 2 — m£) and y = —(A V v^) 2 /m. 

The above model leads to the tri-bimaximal mixing which predicts #13 = 0. It had been the 
focus of A 4 symmetry studies for a few years 13, ll5|, [16]. But it is now ruled out because a non-zero 


has been measured. I 11 this scheme, in order to obtain the tri-bimaximal mixing, the neutrino 
mass matrix with “11” and “33” entries to be equal is crucial. It has been pointed out|13| that a 
more natural form of vev structure will lead to the “33” entry in the neutrino mass matrix to be 
deviate from the “11” entry which leads to a non-zero # 13 . To achieve this, for our purpose here, 
we will introduce two scalars Si> and SV' which are SM singlet but transform as V and 1" under 
A 4 . This results in two new terms for Mr in the Lagrangian 


Ys'{vrVr)i"Si' + Y s »{v r Vr)vSi" + H.C. 


(15) 


After Si'^n develops a non-zero vev, vs>,s ", we have 

( mi 0 m x \ 

0 m 2 0 , (16) 

m x 0 m 3 J 

where m 4 = m + Y S 'Vs r + Y S "Vs", m 2 = m + ou 2 Y S 'Vs' + ouY S "Vs" and m 3 = rn + ojYs>vs> + <jJ 2 Y S "Vs" ■ 
The resulting light neutrino mass matrix M u no longer has w — z, but has 

w = ~Xlv^m 3 /(m 1 m 2 - m \) , z = -X 2 l/ m 1 /(m 1 m 3 - m \) , (17) 

and x and y are changed to 

x = A^m x /(mim 3 - m\) , y = -X 2 v v\lm 2 . (18) 

In the basis where the charged lepton mass matrix is diagonalized, the neutrino mass matrix 
becomes 


1 / w + 2x + y + z w — cu 2 x + co 2 y + uiz w — u>x + ouy + cu 2 z \ 
m u = Uj M U U* — - I w — cu 2 x + ou 2 y + luz w + ‘luux + ouy + ou 2 z w — x + y + z I . (19) 
M w - oux + ouy + ou 2 z w — x + y + z w + 2ou 2 x + ou 2 y + ouz J 

Inserting ou = exp(i2Tr/3) in the above, m u can be transformed into the form in eq.(151) by redehne 
right-handed charged leptons. The parameters in the set Paa : {w,x,y,z} are in general complex 
which will not always have 6 = — 7 t /2 and 0 23 = 7 t/ 4 . One needs to work in the GLS limit which 
can be realized if the parameters in the set Paa are all real. In this case the complexity of the mass 
matrix is purely due to the A 4 group theoretical C-G coefficients ou and ou 2 . This is a case where 
CP violation is caused by C-G coefficients providing a concrete example of intrinsic CP violation. 

Before we analysis the general features of the neutrino mass matrix with complex parameters 
in the set Pa 4 , we would like to analysis the constraints on the model parameters to have the GLS 
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limit, that is, to have w,x,y,z to be real. The complexity of the parameters can appear in the 
Yukawa couplings, in the vevs, and also in places where uu l appear in m,. To make the Yukawa 
couplings and scalar vevs real, one can require the model Lagrangian to satisfy a generalized CP 
symmetry under which 


(*i, il , il) -+ (4) CP , 
($ 1 ,$ 2 ,$ 3 )^($L$S,^) 


(il ) CP , (il) CP ), ( 4 , 

, (Xi , X 2 , Xs) ->• (xl 


4 , 4 ) (( 4 ) CP , ( 4 ) CP , ( 4 ) CP ), 

xS, x!), (20) 


and all other fields transform the same as those under the usual CP symmetry. Here the superscript 
CP in the above indicates that the fields are the usual CP transformed fields. 

The above transformation properties will transform relevant terms into their complex conjugate 
ones. Requiring the Lagrangian to be invariant under the above transformation dictates the Yukawa 
couplings to be real. The same requirement will dictates the scalar potential to forbid spontaneous 
CP violation and vevs to be real. One, however, notices that the parameters m 2j 3 are in general 
complex even if the Yukawa couplings and the vevs of the scalar fields are made real because of the 
appearance of uf. To make them real to reach GLS limit, it is therefore required that 

Irn(uj 2 Ys'Vs' + ojY S "Vs") = Im(ojY S 'Vs' +c o 2 Y S "Vs") = 0 . (21) 

The above can be achieved by the absent of the scalar fields S'" in the theory or set Ys'Vs 1 = 
Ys"Vs"- If the vev structure of x is fixed as given previously, absence of S’ will not have a 
phenomenologically acceptable mass matrix. Therefore, we will take the later possibility as example 
of GLS limit case to show some detailed features. In this case M v can be diagonalized by V u as the 
following 

c 0 — s \ 

010, (22) 
s 0 c J 

where s = sin 6 and c = cos 6. One obtains the mixing matrix to be 

l / c+s 1 c-s \ 

Vpmns = u}v u = —= \ c + us u 2 LUC - s , (23) 

V'J y C + W 2 s CU LU 2 c—S ) 



Normalizing the above mixing matrix to the standard parametrization in eq. (12]h one obtains 


1 1 _ (1 - 2 CS ) 1 / 2 

Sl2 ~ W J S 23 — —7= , Sl3 — -- 

y2(l + cs) v2 v3 

Here we have normalized c i3 and s i3 to be all positive. The neutrino eigen-masses are all real, but 
in general they can take positive or negative values depending on the values of w, x, y and z. Note 
that the absolute values of elements in the second column of Vpmns are all l/\/3. 

We now find the conditions for predicting 5 = —7 t/2 and 6 = +7 t/2 . An easy way of doing this 
is to study the Jarlskog invariant quantity [13] J = ^{VeiV^V^V^)- Eqs.(l2il and (1231) give 



^ ~ c l 3 ' S 12 C 12S23 c 23 s 13 s * n ^ 



(25) 
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which leads to 


-1 , if c 2 > s 2 , 
+1 , if s 2 > c 2 . 


(26) 


S 


71 

= — X 
2 


Note that J is not zero implying CP violation which is caused by the complexity of C-G coeffi¬ 
cients. Eq. (125|) can be transformed into the standard parameterization by multiplying the Vpmns 
on the right and left by diagonal matrices P r = diag( 1, 1 ,i) and Pi = diag( 1, (cu 2 c — s)/\oj 2 c — 
s|, (cue — s)/\ojc— s|), respectively. Pi does not have physical effect because it can be absorbed by 
redefinition of right-handed charged leptons. The physical effects of P r is to change the sign of m 3 . 

Let us now compare experimental data with the model predictions for the mixing angles and CP 
violating phase. There are several global fits of neutrino data[3, HJ. The latest fit gives the central 
values, la errors and the 2 cr ranges as the following!! 


5/IT 

NH l-41±g;|t 

2 cr region 0.0 ~ 2.0 

IH 1.48 ±0.31 

2cr region 0.00 ~ 0.09&0.86 ~ 


'U2 

0.323 ±0.016 
0.292 ~ 0.357 
0.323 ±0.016 
2.0 0.292 ~ 0.357 


4/10- 2 
2.26 ± 0.12 
2.02 ~ 2.50 
2.29 ±0.12 
2.05 ~ 2.52 


^23 

0 567 +0,032 
u.oo /_ 0124 

0.414 ~ 0.623 

D c;7q+0-025 
U.iJ 1 , J_o 039 

0.435 ~ 0.621 


(27) 


Here NH and IH indicate neutrino mass hierarchy patterns of normal hierarchy and inverted 
hierarchy, respectively. In the model above, adjusting the values, w, x, y and z, both NH and IH 
mass patterns can be obtained. There is strong hint that the Dirac phase should be close to 37 r/ 2 (or 
equivalently —tt/2). Therefore one should take the parameter space so that c 2 > s 2 . The value 
—7 t/2 predicted in the model is in agreement with IH within ler range. Although for NH case <5 is 
outside of ler range, there no problem with 2 cr range. For s 2 3 , the model predicts s 23 = 0.5. This 
value is outside of la range for both the NH and IH cases. However, they are, again, in agreement 
with data within 2 cr. 

In the model S 13 = (1 — 2cs) 1//2 /\/3 is not predicted. But one can use information from S 13 
to fix cs = 0.497 ± 0.018 to predict s 2 2 = 0.334 ± 0.004 for both NH and IH cases. This is in 
agreement with data within ler. Note that V 2 2 = (S 12 C 13) 2 = 1/3. It agrees with data within la. 
It is remarkable that neutrino mixing matrix in this model with just one free parameter can be in 
reasonable agreement with data. This may be a hint that it is the form for mixing matrix, at least 
as the lowest order approximation, that a underlying theory is producing. 

If w, x , y and z are allowed to be complex, the GLS is explicitly broken, there are modifications 
to the mixing angles. There is additional source for CP violation other than the intrinsic one 
from complexity of C-G coefficient, and also the mixing angles will be modified. The eigen-masses 
will contain Majorana phases. Detailed analysis of how to diagonalize the mass matrix has been 
discussed in Ref.jl2[. In general this model does not always predicts 5 = ±7t/2 and 623 = 7t/4. The 
mixing matrix can be, in general, written as 


Vpmns 


U\V P K = 


1 

7! 


c ± se lp 1 
c ± c ose ip cu 2 
c ± cu 2 se* p u 


ce lp — s \ 
Luce lp — s 

l o 2 ce lp — s J 


(28) 


where V p is a diagonal matrix diag( 1,1, e ip ) with tanp = Im(xw* ± x*z)/Re(xw* ± x*z). It is 
interesting that the phase p does not show up in J which is still — (c 2 — s 2 )/ 6\/3- This implies that 
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CP violation related to neutrino oscillation is still purely due to intrinsic CP violation. The mixing 
angles and the Dirac phase 5 are all modified with 


Sl2 = 


\/ 2 (l + cs cos p ) 1 / 2 


S23 — 


(1 + cs cos p + \/3 cs sin p ) 1 4 / 2 
a/2(1 + cs cos p )5 


. u = d— 2 c 7^> 1/2 . (29) 
v 3 


and 


r 4c 2 s 2 sin 2 p. 1/9 . 3c 2 s 2 sin 2 p . 

Sln<5 = t 1 + 73-^5 t)“ ' t 1 - TTVWvAj) 


1/2 


(1 + cs cos p ) 5 


X 


-1 , if c 2 > s 2 , 
+1 , if s 2 > c 2 . 


(30) 


In this case, the new parameter p can be used to improve agreement of the model with data. In 
both NH and IH cases, <5 and S 23 can be brought into agreement with data at lcr level. To see how 
this can be done, as an example, we take the largest value of cs so that S 13 takes its lower la allowed 
value, and then varying cos p to obtain the upper la allowed value. This fixes cs and cos p to be 
0.468 and 0.992, respectively. With these values, S 23 and 6 are determined to: 0.534 and 1.4267T, 
respectively. These values are in agreement with data at la level. When more precise experimental 
data become available, the model with complex model parameters can be distinguished from that 
with the parameters are all real and other models. 

In summary we have shown that neutrino mass matrix reconstructed with 5 = —7t/ 2 and 
$23 = tt/ 4 has several interesting properties. We find that a theoretical model based on the A 4 
symmetry naturally realize the GLS limit and predicts such a neutrino mixing pattern together 
with the prediction \V e 2 \ = l/\/3. In this model, CP violation can be solely come from the complex 
group theoretical C-G coefficients if the neutrino Majorana phases are zero or n. This model fits 
experimental data very well and can be taken as the lowest order neutrino mass matrix for future 
theoretical model buildings. If there are additional source of CP violation other than those intrinsi¬ 
cally existed in the C-G coefficients, the CP violating phase 6 and the mixing angle #23 can be away 
from —7 t/ 2 and 7 t/ 4. The models discussed can fit data within la. Future improved experimental 
data will be able to further test the model and provide more hints for the underlying theory of 
neutrino mixing. 
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